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Abstract 

Interactions of ultrahigh energy neutrinos of cosmological origin in large 
volumes of dense, radio-transparent media can be detected via coherent 
Cherenkov emission from accompanying electromagnetic showers. Antarc¬ 
tic ice meets the requirements for an efficient detection medium for a radio 
frequency neutrino telescope. We carefully estimate the sensitivity of realis¬ 
tic antennas embedded deep in the ice to 100 MHz - 1 GHz signals generated 
by predicted neutrino fluxes from active galactic nuclei. Our main conclu¬ 
sion is that a single radio receiver can probe a ~ 1 km^ volume for events 
with primary energy near 2 PeV and that the total number of events regis¬ 
tered would be roughly 200 to 400 year“^ in our most conservative estimate. 
An array of such receivers would increase sensitivity dramatically. A radio 
neutrino telescope could directly observe and test our understanding of the 
most powerful particle accelerators in the universe, simultaneously testing 
the standard theory of particle physics at unprecedented energies. 



1 Introduction 


An ultrahigh energy (UHE) neutrino telescope is a multipurpose instrument 
by its very nature. It has the potential to probe, in a fundamentally new 
manner, the most distant and energetic objects in the universe (such as 
active galactic nuclei or AGN) as well as exotic structures within our own 
galaxy [1], Due to the enormous energies of the particle interactions involved 
(> lO^^eV), it also serves as a natural laboratory to test our understanding 
of fundamental physics in the UHE regime. 

As a new observational window on the universe, neutrino astronomy’s ul¬ 
timate impact on our understanding of astrophysical phenomena cannot be 
known in advance. The history of astronomical science is certainly encourag¬ 
ing: every new observational tool has produced unforeseen discoveries. In the 
case of neutrino astronomy, the prospects seem to be particularly promising 
since for the hrst time one would use a fundamentally different particle to 
view the cosmos. Unlike photons, which are easily stopped by small amounts 
of intervening matter, neutrinos interact very weakly, as if custom-made to 
give us information about regions of the universe which can not be directly 
probed using the techniques of conventional astronomy. Of course, this same 
property also makes neutrino observation a technical and scientihc challenge 
[ 1 , 2 ]. 

As an example of the particle physics aspect of a neutrino observatory, 
consider the fundamental neutrino-nucleon cross section, . When the 
invariant center of mass energy, is much larger than the W-boson mass, Mw, 
the magnitude of this cross section is directly determined by the number of 
quarks carrying a small fraction, x, of the total nucleon momentum in the 
center of momentum frame [3, 4, 5, 6]. The quark distribution comes from 
fundamental, non-perturbative, strong interaction processes. It follows that 
UHE neutrino detection probes not just the depths of the cosmos but also 
intricate details of hadronic structure and strong interaction physics. One 
can only speculate about what other interesting phenomena might present 
themselves at untested energy extremes. Eor example, it has been suggested 
[7] that tan neutrinos resulting from neutrino oscillations would produce 
a clear signal in a UHE detector, providing evidence for “new physics” if 
observed. 

Several new high energy detectors utilizing optical sensor technology are 
now coming into operation. These operate mainly by detecting muons from 
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charged-current muon-neutrino interactions. The muon’s path must inter¬ 
sect a detection volume arrayed with photomultiplier tubes to amplify visible 
Cherenkov photons the muon emits. The long range of the muon is helpful 
in probing a large target volume. This detection strategy appears to be ideal 
for observing neutrinos in the GeV to TeV range and it is being employed 
in a variety of naturally occurring, dense and optically transparent environ¬ 
ments such as deep freshwater lakes (NESTOR, Lake Baikal), deep ocean 
(DEMAND), and Antarctic ice (AMANDA) [1]. However, the combination 
of flux and detector efficiency severely limits this approach in the energy 
regime above 100 TeV. 

Here we consider an attractive alternate detection strategy based on 
coherent radio Cherenkov emission from neutrino-induced electromagnetic 
showers [8, 9]. This approach was hrst developed and tested in the context 
of cosmic ray induced air showers [8, 9, 10]. The advantages of going to 
a dense radio-transparent medium are spectacular but have only been fully 
appreciated relatively recently. The radio power emitted coherently scales 
inversely with the radiation length of the medium squared. About 1 million 
times more radio energy is emitted by a shower in ice than in the old air 
shower experiments. Cold ice happens to be an abundant natural substance 
that is quite transparent to radio signals, having an attenuation length of 1 
km or more at Antarctic temperatures. Thus, Markov and Zheleznykh [11] 
suggested using the Antarctic ice as a solid natural target for neutrino detec¬ 
tion, and this idea has been subsequently developed by a number of authors 
[12, 13, 14, 15]. 

The radio power emitted coherently also goes like the primary shower 
energy squared [8, 9]. (This can be compared with emission from charged 
shower particles that add incoherently as in the case of visible light, which 
rises at most linearly with energy.) The power received at a distance r from an 
event falls like 1/r^ from geometry (assuming attenuation can be neglected). 
As a rule of thumb, for radio the detection range increases proportionally to 
the primary energy. This is a very interesting effect. A range increasing lin¬ 
early with energy translates into a target volume going like the energy-cubed. 
Earlier numerical estimates including realistic antennas and signal/noise = 
1 found radio signals detectable at a distance as large as 1 kilometer if the 
energy were 1 PeV [12, 15]. By the scaling laws, one could also detect 100 
TeV neutrinos as far away as 100 m, or 10 TeV neutrinos at 10 m. It follows 
that radio detection is not limited intrinsically to fantastically high energy 
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thresholds. What radio does is offer a way to probe enormous target vol¬ 
umes if UHE neutrinos exist. For PeV neutrinos and above, radio can probe 
cubic kilometer-scale volumes per detector element in the near future at a 
comparatively low cost. 

Considerable attention is being given now to “KM3” proposals for visible- 
light-based detection with target volumes of the order of Ikm^ [1, 16]. These 
proposals are ambitious and exciting. Due to the different scaling laws with 
energy the radio and optical approaches are complementary rather than in 
direct competition. Depending on the primary flux, we And that the energies 
at which radio detection is most exciting range from 100 TeV to 100 PeV, 
bridging a gap in conventional methods and considerably extending the range 
of energy. It is a potentially valuable coincidence that the medium chosen for 
the AMANDA optical array (pure —50°C ice) is also beautifully transparent 
to radio signals up to the GHz range [17] where coherence is maintained. 
The benefits of integrating radio and optical detection seem to be very great 
indeed and should be carefully considered. 

It might seem that a few radio antennas placed on the ice surface, either 
alone or in conjunction with a conventional array, would be a good way 
to extend existing studies to above the PeV range. Early calculations by 
Ralston and McKay [13] assumed this geometry and detection of upward- 
going showers ignoring attenuation. However, Earth-shadowing effects cut 
the flux of upward going neutrino events above 100 TeV severely [18]. A 
reliable estimate must include the angular dependence of the attenuation 
[15]. The flux of of “horizontal” neutrinos which penetrate with minimal 
attenuation actually dominates. We present here the first calculations in 
which the full effects of the geometry are taken into account and a variable 
antenna depth is included; we consider optimizing the detection by studying 
buried detectors with good horizontal acceptance as a function of depth, to 
take advantage of unattenuated sideways flux at 1 PeV and above. Buried 
detection has not been considered in detail before. The results represent a 
substantial improvement in the whole scheme. 

It is natural to contemplate a single telescope facility employing both op¬ 
tical and radio detection to observe cosmic neutrinos over an enormous range 
of energies. However, radio detection has enough advantages to make an inde¬ 
pendent facility very attractive. We And several hundred events from electron 
neutrinos, per detector per year, using standard flux estimates. Backgrounds 
from atmospheric neutrinos are negligible. These topics are discussed in de- 
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tail in the sections which follow. 

Because of the unusual growth in detection efficiency with energy, we have 
focused on the energy region of 100 TeV and above. In this regime, much of 
the recent activity in neutrino astronomy has been stimulated by the prospect 
that AGN may be intense sources of energetic neutrinos [19, 20, 21, 22]. 
Theoretical results for “typical” AGN have been integrated to estimate the 
diffuse neutrino flux to be expected from all the distant active galaxies in 
the universe. In this work we estimate the sensitivity of radio detection in 
Antarctic ice to two different predictions for this isotropic or unresolved AGN 
flux. Our results show that even a single radio receiver is capable of detecting 
AGN at the current predicted flux levels. 

Not only would AGN be detected, but we might also be able to distin¬ 
guish between competing models of the underlying physics because of the 
different energy dependences for the predicted flux. Very generally, if a pri¬ 
mary neutrino integral spectrum is falling as with an unknown spectral 
index 7 , and the efficiency for detection V(E) is a known, increasing func¬ 
tion of energy, as it is in our case, then the distribution of events, , 
will have a characteristic maximum at some energy Amax- This permits a 
direct measurement of the unknown source spectrum. We will show that the 
unknown spectral index can be found from an observed value of Amax by 
evaluating 7 = Umax- the reader will see in subsequent sections, the 

two AGN models we consider have spectral indices that differ by roughly 1 
[19, 20, 21]. Given our evaluation of the efficiency for radio detection, we 
predict an order of magnitude difference in the value of Amax between the 
two models. With even modest energy resolution, one could reasonably hope 
to distinguish between these two pictures of AGN physics. 

Figure (1) is a sketch indicating the geometry of a typical detected neu¬ 
trino event to help the reader follow the sequence of concepts and calculations 
developed in the various sections of this paper. Sections 2 and 3 present esti¬ 
mates for the UHE flux incident at the Earth’s surface and a detailed account 
of how this nominal flux is modihed by interactions as it passes through the 
Earth on its way to the region containing the detector. The sketch indicates 
a typical UHE neutrino that enters the detection region and has a charged- 
current interaction with a target nucleon; a calculation of the rate for these 
interactions is the topic of Section 4. The resulting charged hnal state lepton 
initiates an electromagnetic shower of charged particles and photons. The 
production rate of these cascades is the subject of Section 5. Electromag- 
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Figure 1: The geometry of a typical detected neutrino event. Our calculation 
includes all known effects: particle physics in the fundamental interaction 
cross sections, the angular dependence of attenuation in the Earth, electro¬ 
magnetic shower evolution, coherent Cherenkov radio emission, polarization 
of the emitted held, pulse dispersion in the ice and antenna characteristics. 


netic cascade evolution produces a net excess of electrons compared with 
positrons that gives a coherent Cherenkov pulse at radio frequencies. This 
pulse is a ring-shaped structure, indicated in the sketch by dashes that prop¬ 
agates along the surface of a cone dehned by the Cherenkov angle for the 
medium. The characteristics of the radio signal and a detailed treatment of 
its propagation from the site of the electromagnetic shower to a point where 
it intersects the antenna array is the topic of Section 6. In Sections 7 and 
8 we present an analysis of how the radio pulse couples to a realistic an¬ 
tenna/receiver combination and how it is detected against a model for the 
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noise background. Using a signal-to-noise criterion, we present in Section 9 
a calculation of the effective detection volume per antenna that incorporates 
all the threshold effects relevant to detection. The calculations of shower 
rates and effective volumes come together in Section 10 where we show the 
detected number of events per year per antenna. In Section 11 we give a 
procedure for determining the source spectral index from experimental data. 
Finally Section 12 summarizes our results and estimates rates for a multi-unit 
neutrino telescope. 


2 Predictions for Diffuse UHE Neutrino Fluxes 

Figure (2) shows recent predictions for the unresolved AGN neutrino flux 
due to Szabo and Protheroe (denoted hereafter as SP) [20] and Stecker et 
al. (denoted ST) [19]. The Szabo and Protheroe flux we show is the most 
optimistic of their range of estimates. We have chosen these two AGN flux 
models as representative of the range of current theoretical prediction. Also 
shown is the atmospheric neutrino spectrum as predicted by Lipari (denoted 
AT) [23]. The the curves are piecewise linear approximations (on a log-log 
plot) to the published curves used for convenience. These three spectra will 
serve as the flux input for our calculations. The hgure depicts the -|- 77^ 
flux at the Farth’s surface. The ratio of electron neutrinos to muon neutrinos 
is approximately 2/3 in the SP model and 1/2 in the ST case. 

We adopt the following notation. An isotropic differential neutrino flux 
at the Farth’s surface can be expressed, 

^{E.) = UEo) * • ( 1 ) 

Between the Farth’s surface and the detector, these neutrinos must pass 
through different amounts of intervening matter. Cross sections for charged 
and neutral current interactions are sufficiently large at the energies of in¬ 
terest that we must ask how the originally isotropic flux appears after its 
journey. 
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Figure 2: Three neutrino flux models used in the present calculati 
spectra for diffuse AGN are due to Szabo and Protheroe (SP) an( 
et al. (ST). The third curve is the atmospheric neutrino prediction 
lATi. 




3 Modification of Flux by Its Passage Through 
Earth 

The neutrino flux which reaches a buried detector is modified by interactions 
with matter as it passes though the Earth. Aside from the reaction, I7g + e“ —> 
W~ —> anything, which is large only near the resonant energy 6.4 PeV, 
scattering of neutrinos on atomic electrons can safely be neglected, leaving 
charged and neutral current interactions with nucleons as the two relevant 
mechanisms to consider. The flux is modifled in two distinct ways. Charged 
currents, ly -h p —A" + e“ and 77 + n —X e~^ directly remove neutrinos 
from the flux while neutral currents, u p —X u and 77 + n — A" + 77 
shift neutrinos into lower energy bins. 

The effect is energy dependent owing to the energy dependence of the 
relevant interaction cross sections. We have recently calculated the charged 
current neutrino-nucleon cross section at ultra-high energy [ 6 ]. Our calcu¬ 
lation is based on QCD evolution of recent electroproduction data at small 
Bjorken x from the HI and ZEUS collaborations at HERA [24, 25]. Eig- 
ure (3) represents our current knowledge of the cross section over the GeV 
to PeV range showing previous estimates [3, 4, 5] for the UHE region that 
were based on incomplete knowledge of the small x physics along with our 
improved calculation [ 6 ]. The region of existing high energy cross section 
data is indicated, along with one very high energy data point near 50 TeV 
extracted by the HI group who measured (ri^ep —> u-\- X) [26]. The lone data 
point is at approximately 10 times the highest energy previously reported 
and is in good agreement with our new result, thus tying together two com¬ 
pletely independent aspects of HERA measurements. Our result is roughly 
2.2 times previous estimates for neutrino energies of 1 PeV and above. The 
following parameterization of our result is valid at the 5% level of accuracy 
over the range 50 TeV < E < 50 PeV, 

= (1.03xl0-3^m2)(;|^)°-«22exp[-0.0231 ln^(;|^)] • (2) 

At these large energies, the cross section is dominated by events in which 
the struck parton carries a very small fraction of the total nucleon momentum 
(he. small x). In this regime, the quark sea dominates over valence contri¬ 
butions and the equal abundance of quarks and anti-quarks means that the 


8 



Cross Section (cm^2) 
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Figure 3: Total charged current neutrino-nucleon cross section versus incident 
neutrino energy. The extrapolated 4-Fermi result is shown along with the 
single high energy HFRA data point. Also shown is the new ultrahigh energy 
prediction of Frichter, McKay and Ralston [ref. 6] along with previous results 
of McKay and Ralston [3] and Reno and Quigg [4]. 
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neutral and charged current cross section are proportional. The constant of 
proportionality is related to the weak mixing angle, a well known quantity. 

^ O.W afciE.) (3) 

In addition to the energy dependence, the modihed flux gets an angular 
dependence due to differing amounts of matter it must traverse en route to 
the detector. The overall effect is expressed by an evolution equation for the 
flux. Consider charged and neutral current events occuring over an interval, 
dt = n]\j{z)dz, with n]\j{z) being a local nucleon density for some position 2 ; 
along the path. 


dXncf) 

dt 


(E.) 






dE', 


(t){E^) dE^ 


{El,E, 


(4) 


The effective neutrino-nucleon cross section, is a sum of charged and 

neutral current contributions evaluated at Cj,, minus a correction which ef¬ 
fectively reduces the cross section at the energy of interest because of neutral 
current interactions which occur at larger energies and have a hnal state 
neutrino energy Eu. Since the magnitude of the ‘T regeneration” correction 
depends on the number of neutrinos available at energies larger than ifj,, the 
effective cross section will be flux dependent via the spectral shape, ^f{Eu). 

Now let 9 be the nadir angle (the angle measured from downward vertical) 
of a particular primary neutrino momentum, and n]\j{z,9) be the nucleon 
density along the path as a function of 9. The effects of the evolution can 
then be expressed simply in terms of the effective cross section [27], 


(t>{E^,9) = (/)o(£'^)exp[-crefr(£'^,7)t(6')] 

= UE.)KE..0), (5) 

with t{9) = / n]\j{z,9)dz and here (j)o{E^) is the incident neutrino flux of 
Eq.(l). We make use of the PREM Earth structure model [28] to obtain 

m- 

In Eig.(4) we show effective cross sections for the SP and ST flux predic¬ 
tions along with the charged plus neutral current cross section to show the 
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effect of the v regeneration term of Eq.(4). The effect is largest for the ST 
flux model at energies below roughly 1 PeV owing to the spectrum’s flatness 
compared with the SP flux. Higher energy neutrinos feed the lower energy 
flux more efficiently in the ST case than the SP case. 

With the effective cross sections in hand, we can use Eq.(5) to hud the 
angular variation of the flux in the region of our detector. Eigure (5) shows 
our result for the SP and ST AGN flux models. The angular variation of the 
attenuated flux clearly demonstrates the effect of shadowing by the Earth’s 
mass which becomes ever more pronounced with increasing energy due to 
the rising interaction cross section of Eq.(2). In the calculations which fol¬ 
low, we assume that the flux as given by Eq.(5) is invariant with respect 
to translations within several kilometers of the detector. This is reasonable 
since such distances are very small compared with a typical mean free path 
for the neutrinos. 


4 Charged-Current Lepton Production In De¬ 
tection Region 


Using the local neutrino flux, we can directly evaluate the number of leptons 
produced by charged-current interactions in the vicinity of the detector. We 
let (/)® and represent the local neutrino fluxes of electron and muon type as 
determined by Eq.(5). The rates, Pe,^, for electron and muon production in 
the neighborhood of the detector are found by integrating over all potential 
parent neutrinos having energy, Uj,, weighted by the appropriate differential 
cross section for production of charged leptons with energy, Ee or E^\ 


1‘co 


( 6 ) 


Our results are shown in Pig. (6) and indicate the number of leptons produced 
per year, per energy bin, per unit of solid angle, per unit volume. Under the 
assumptions of the SP flux model, 40% of these are and the remaining 
60% are /i^. In the ST case, the corresponding percentages are 33% and 66%. 
Each charged lepton has the potential to initiate an electromagnetic (EM) 
shower in the ice that is detectable by virtue of its accompanying Cherenkov 
shock wave. Our evaluation of the EM shower rate is the topic of the next 
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Nadir Angle (deg) Nadir Angle (deg) 

Figure 5: Diffuse AGN neutrino flux in the detection region showing calcu¬ 
lated angular variation due to Earth shadowing. Neutrino flux versus angle 
with respect to vertical downward is plotted for a series of neutrino energies 
from 100 TeV to 10 PeV. Curves are labeled with log^Q[energy(TeV)]. Parts 
(a) and (b) show the total flux of i/g, G, and using the SP and ST 
predictions respectively. 
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(b) ST 




Nadir Angle (deg) Nadir Angle (deg) 

Figure 6: Charged lepton production per year per TeV per steradian 

per cubic kilometer of polar ice due to charged current AGN neutrino interac¬ 
tions. The production rate versus angle with respect to vertical downward is 
plotted for a series of neutrino energies from 100 TeV to 10 PeV. Curves are 
labeled with log^Q[energy(TeV)]. Parts (a) and (b) show the rates induced 
by the SP and ST flux predictions respectively. 
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section. 


5 Electromagnetic Shower Rates 

We will present detailed calculations for radio detection of the process Ue + 
p —> e +W, followed by an electromagnetic shower caused by the energetic h- 
nal state electron. We ignore radio energy produced by the hadronic shower, 
whose evolution involves strong interactions. It follows that our calcula¬ 
tions of radio power and event rates are very conservative. Provorov and 
Zheleznykh [15], for example, treat hadron showers on the same footing as 
electron showers. The radio power from Ue induced showers is nearly doubled 
by this effect, and events produce a full electron-like shower. Additionally, 
neutral-current neutrino interactions also produce observable showers. We 
will return to this point in the hnal section of this paper. 

Contrary to optical detection, we hud that electrons are more efficiently 
detected with radio than muons. This fact holds despite the long range of 
the muon. Every electron produces a shower carrying all of its energy, and 
whose compact evolution produces a sizable radio signal. Thus we have for 
electron showers the simple relation 0) = Ve{Ee-,G)- On the other 

hand, it is rare for a muon to bremsstrahlung a photon carrying a sizable 
fraction of its energy which would make for efficient radio detection. The 
small probability for bremsstrahlung closely compensates the effects of the 
muon’s long range. Moreover, for every photon carrying a fraction x of the 
muon’s energy that will produce a shower carrying energy Eg, the muon has 
to have energy E^ = Egjx. The distribution over x tends to force E^ up into 
a region of higher energy and lower flux, decreasing the rate. 

To be more quantitative, we adopt a continuous loss expression for the 
muon shower rate 0) for showers moving at angle 6 with respect 

to the nadir, given by 

rco Pj^ rX[0) 

rp--{Eg, e) = Eg) e) , (7) 

with 

K = (E. + - £ . ( 8 ) 

where e and ^ are the usual energy loss parameters, and is the bremsstrahlung 
probability [29]. The range integration of Eq.(7) hnds the muon flux at en- 
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ergy arising from muons produced at range r and having an initial energy 

given by the range-energy relation Eq.(8). The upper integration limit, 
X{0), is the chord length from the detector site to the Earth’s surface for a 
hxed direction 9. 

Eigure (7) shows our result for the number of EM showers occuring per 
year, per energy bin, per unit solid angle, per unit volume due to muon 
bremsstrahlung events. The curves have an interesting shape which peaks in 
the horizontal direction. This is to be expected since upward going muons 
are few in number due to neutrino absorption by the Earth, and smaller 
target volumes above the detector produce fewer muons compared with the 
more favorable horizontal direction. Recalling from Section 2 that muon 
neutrinos outnumber electron neutrinos by 3:2 in the SP case and 2:1 in 
the ST case, we nevertheless hud that the peak horizontal shower rate is a 
bit more than a factor of 10 lower than the rate induced by the electron 
neutrino interactions. Given that the contribution of muons is small, the 
continuous energy loss expression, which neglects stochastic fluctuations of 
photon emission, is quite adequate for our study. 

We now want to determine the number of these showers giving rise to 
a detectable signal in a radio receiver located some distance below the ice 
surface. This task involves several steps that will be dealt with in subse¬ 
quent sections. Our calculation takes into account the details of EM shower 
structure in ice, a characterization of the resulting radio pulse, pulse attenu¬ 
ation as it propagates from interaction site to antenna, the coupling of signal 
to antenna voltage including polarization effects, the receiver characteristics 
transforming the induced antenna voltage into a detectable signal, and hnally 
an estimate of the noise levels with which the desired signal from neutrino 
interactions must compete. 

In the sections which follow, we will fold together what is known about 
each of these items and arrive at a conservative estimate for detectible neu¬ 
trino events. This will then be expressed in terms of an energy/direction- 
dependent (E, 9) sensitive volume for an individual sensor (antenna plus 
receiver). This effective volume can be directly multiplied by and 

p^hower arrive at an event rate predictions for electron and muon neutrinos 
respectively. In this way threshold effects are conveniently incorporated into 
the overall rate estimate. 
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Nadir Angle (deg) Nadir Angle (deg) 

Figure 7: Number of EM showers induced by muon bremsstrahlung events 
per year per TeV per steradian per cubic kilometer of polar ice. The pro¬ 
duction rate versus angle with respect to vertical downward is plotted for a 
series of neutrino energies from 100 TeV to 10 PeV. Curves are labeled with 
logio[energy(TeV)]. Parts (a) and (b) show the rates induced by the SP and 
ST flux predictions respectively. 
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6 The Cherenkov Pulse and Its Propagation 
in Antarctic Ice 

Scaling laws governing the overall behavior of radio detection have been un¬ 
derstood for some time now [9] based on analytic models of shower structure 
together with the characteristics of Cherenkov emission. Recently, a signih- 
cant advance in our understanding of coherent radio pulses from showers in 
ice was due to the detailed numerical study of Zas, Halzen and Stanev [14], 
Their Monte Carlo code incorporates all the relevant particle interactions, 
tracks the history of each shower particle and then coherently sums the low 
frequency Cherenkov radiation from all the resulting charged particle tracks. 
The simulations give a very detailed view of the shower as it propagates 
through the ice with its accompanying radio pulse. 

The relevant result for our purposes is the variation of the signal frequency 
spectrum with primary shower energy and also the angle between the longi¬ 
tudinal shower axis and the observation direction. Near the Cherenkov angle 
{9c = 56° in ice), the Fourier transform of the pulse electric held, E{u^ i?, 9c), 
can be parameterized compactly as follows [14], 


R\E{v, R,9c)\ 


0.55 X 10 '^{uluo) E,s 1 9 - 9c ly . 21 / ^ 

1 + 0.4(i//i/„)2 lTeV®^^^~2^ 9c, ^^MHz 


(9) 


where R is the distance between observer and cascade, Uo = 500 MHz, 
9o ~ 2.4°, and Eg is the total shower energy. The deliverable power in 
this signal is proportional to the square of the electric held and scales as 
the primary energy squared with a peak near 1 GHz. This E^ growth of 
the available signal power is to be contrasted with an Eg dependence for 
incoherent emission at higher frequencies. 

The next step is to determine how interaction with the ice modihes the 
pulse dehned in Eq.(9) as it propagates from event site to antenna. Below 
about 5 GHz, dispersive ehects on the pulse can be neglected, leaving ab¬ 
sorption of pulse energy as the relevant ehect to consider. The absorption 
is frequency and temperature dependent. We use the ice temperature versus 
depth data taken at the Antarctic station Vostok by Salamantin et al. [30] 
along with measurements of attenuation length variations with frequency 
and temperature due to Bogorodsky and Gavrilo [17] to predict the pulse 
spectrum after propagation along a given path from event to receiver. 
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Figure 8: Expected absorption length in kilometers versus depth and fre¬ 
quency based on Antarctic temperature measurements at the Vostok station 
[ref. 30]. Contours are in increments of 0.2 km from 1 to 4 km. 
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In Figure (8) we show the resulting attenuation length as a function 
of frequency and depth below the ice surface. At 400 meters depth, the 
absorption length ranges from 4.4 km at 100 MHz to 2.0 km at 1 GHz. 
Just above 1 GHz new loss mechanisms become active and the attenuation 
increases sharply. Because it has the lowest temperature, the ice nearest 
the surface is the most radio-transparent. By itself, this factor would suggest 
shallow deployment of a radio receiver for optimum event detection. However, 
recall from the results of Section 3 that most of the neutrino flux is coming 
from the upper hemisphere due to the phenomenon of shadowing by the 
Earth. This means that there are advantages to having a large detection 
volume above the receiver: that is, deep deployment. The depth which yields 
the optimum overall rate of detectable events is a quantity that is vital for 
detector design, and we will evaluate it in subsequent sections. 

As an example of signal attenuation. Fig. (9) shows attenuated spectra for 
pulses propagating a distance, R = 2 km, to a receiver located at a depth 
of 700 m from a series of hve angles ranging from 30° to 110° measured 
with respect to vertical downward. The variation of the spectrum magnitude 
and corresponding pulse height with angle underscores the virtues of lower 
temperature layers and the advantages of operating the detector in the more 
radio-transparent ice. 

7 Coupling of Signal to Detector 

The attenuated pulse from the previous section arrives at the antenna having 
a particular polarization determined by the geometry of the event. The elec¬ 
tric held, A, is transverse to the propagation direction and lies in the plane 
containing the pulse’s Poynting vector and the longitudinal shower axis. This 
held induces an open-circuit voltage, 14, in the receiving antenna that can be 
expressed in terms of an ehective height vector, h = h(/)4, which, together 
with antenna gain, expresses the efhciency of the antenna as a transducer 
between the propagation medium and a transmission line connected to the 
antenna terminals. The angular variable, 4, refers here to the angle between 
the incoming signal’s Poynting vector and the symmetry axis of the antenna. 
The ehective height function depends on the material properties and geome¬ 
try of the antenna itself as well as the dielectric properties of its environment, 
in this case the ice in which it is imbedded. The relation is 
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Figure 9: Cherenkov pulse after propagating 2000 meters from the cascade 
site to a receiver at 700 meters depth. The pulse is shown in the frequency 
domain (modulus only) and time domain with the curve labels indicating the 
nadir angle in degrees for the event. Pulses passing through the colder ice 
layers suffer less attenuation. 
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Figure 10: An equivalent circuit representing the antenna-receiver combi¬ 
nation. These quantities will be used in the evaluation of noise and signal 
charateristics at the receiver output. 


v:= ^€E(f.i’)-h(f)G{i’). (10) 

where G{'tp) is the gain of an omnidirectional antenna, and = 1.78 is the 
refractive index of ice. 

The voltage of Eq.(lO) drives the equivalent receiver plus antenna cir¬ 
cuit depicted in Fig.(10), where Za{f) and ZL{f) represent the frequency 
dependent complex impedances of the antenna and the receiver load re¬ 
spectively. Under matched impedance conditions, the receiver is such that 
ZL{f) = Z*{f) ensuring maximum power transfer from antenna to receiver. 
In this case, the signal voltage we measure across the receiver load is simply. 
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For our numerical work we will use theoretical predictions for the effec¬ 
tive height, gain and impedance of a particular wide-angle biconical antenna. 
Expressions for these antenna characteristics can be found in references [31]. 
Existing impedance measurements for conical antennas [32] match the theo¬ 
retical values quite well. We model a biconical antenna having a total length 
of 22 cm and a conical half-angle of 30°. Results are shown in Fig.(11). One 
can see from the effective height that the sensitivity peaks in a band from 
roughly 100 to 400 MHz and falls off sharply at higher frequencies. At hrst 
this lack of efficiency near 1 GHz might seem like a serious liability since the 
initial pulse spectrum of Eq.(9) peaks there. Recall, however, that the signal 
spectra shown in Fig. (9) peak at lower and lower frequency as attenuation 
is increased. This means that signals near the detection threshold are likely 
to have their power concentrated at frequencies where our biconical anten¬ 
nas are most sensitive. Most of our detectable events will come from large 
distances (an volume effect) and near threshold. From the point of view 
of maximizing the rate of detectable events, the sensitivity of the biconical 
described above seems to make it a good choice. 

The biconical antenna should perform somewhat better than, for exam¬ 
ple, a simple dipole of the same length because of its larger bandwith, helping 
to discriminate the signal pulse from noise. R is also an attractive option be¬ 
cause it has a maximum horizontal diameter of about 10 inches, a size already 
shown to be practical for placing down deep bore holes by the AMANDA 
collaboration. 

The impedance and effective height functions of Fig.(11) are used in Eqs. 
(10) and (11) along with the attenuated pulse spectrum from Section 6 to 
determine the signal at the receiver output. As an example. Fig.(12) shows 
the spectrum at the receiver output, 14(/), for the hve signals described in 
Fig.(9). This is the signal we wish to detect. 


8 Signal to Noise Ratio 

The total voltage observed at the receiver’s output can be expressed as the 
sum of the signal given by Eq.(ll) and a random thermal noise component. 
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Figure 11: Antenna driving point impedance and effective height versus fre¬ 
quency. These are theoretical predictions for a 30° (cone half-angle) biconical 
having a 22 cm total length. 
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This includes Johnson noise generated in the circuitry of the receiver itself, 
as well as other random electromagetic holds present in the environment that 
are absorbed by the antenna. This type of noise is characterized by a hat 
spectrum whose magnitude can be expressed in terms of an equivalent noise 
temperature, r„. The noise power in a band. A/, is simply kTnAf. We 
assume that detection of the neutrino signal will be limited by thermal noise 
and that backgrounds such as radio emission from the galactic plane or man¬ 
made sources are not important. This assumption seems reasonable based 
on measurements of ambient radio noise in the Antarctic by Boldyrev et al. 
[33]. These measurements indicate that the value r„ = 300K is a reasonable 
estimate over the 100 MHz to 1 GHz band that we will be concerned with in 
this work. 

Detection of a known signal in a white Gaussian noise background is 
a standard problem of signal processing and is optimized in the sense of 
maximizing the signal to noise ratio, 

SNR output signal power 

output noise power 

by use of a matched correlation receiver. The matched hlter employs a re¬ 
ceiver transfer function that is ‘matched’ to the complex conjugate of the ex¬ 
pected receiver input voltage. In principle we know the form of the expected 
signal and can make use of the increased discrimination against thermal noise 
that this procedure permits. Under matched conditions, the signal to noise 
ratio can be expressed, 

^ 2 (/-"g' mnMr 

where Rl is the resistive component of the load impedance, Z^. A discussion 
of thermal noise, the matched hlter and the associated expression for SNR 
can be found in reference [34]. We will consider detectable any event having 
SNR > 1. This is a very conservative method because use of correlations 
between different detectors can improve upon this considerably. 
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9 Effective Volume for the Matched Bi-Conical 
Antenna 

Consider a fixed shower energy, Eg, incident neutrino direction, 0, and a 
variable position for the electromagnetic shower relative to the antenna. The 
collection of potential shower positions that satisfys SNR > 1 using Eq.(13) 
dehnes an effective volume for events with hxed [Eg, 9). This volume can 
then be multiplied by or to directly give the rate for events 

initiated by neutrinos of electron and muon type respectively. 

As an example. Fig.(13) depicts a vertical cross section of ice containing a 
receiver 900 meters deep. The contours indicate SNR=1 for shower energies 
of Eg =1 PeV and 100 PeV with the incoming neutrino direction, 6 = 90°, 
hxed (horizontal hux coming from the right side of the hgure). At each en¬ 
ergy, the two lobes (the upper lobe is cut off by the surface) are sections of the 
complete effective volume that can be imagined by rotating the lobes about 
an axis containing the receiver and pointing in the direction of the incident 
neutrino hux. This volume is simply an ‘image’ of the signal Cherenkov cone 
(with some distortions introduced by the variable properties of the ice with 
depth), truncated by its intersection with the ice surface at energies above 
about 1 PeV. This “fat” cone-shaped region is the ehective volume of the 
receiver for that particular energy and hux direction. Each possible combi¬ 
nation of Eg and 6 has its own cone-shaped volume. Note that at energies 
above 1 PeV the lobes extend more than a kilometer from the receiver. One 
could make similar sets of contours for other values of 6 and receiver depth 
in order to see what sections of the ice are probed by a particular placement 
of the antenna. This procedure would be useful in designing a radio receiver 
array. 

Figure (14) shows results for the ehective volume of a 600 meter deep 
antenna versus nadir angle (the angle with respect to local downward vertical) 
and total cascade energy. Curves are drawn for a series of energies from 0.25 
to 32 PeV. Considering the energy dependence of Eq.(9) along with the 
conical geometry of the sensitive region, one predicts a volume scaling like 
E^ in the absence of attenuation. In practice, both the signal loss and the 
cutoh presented by the ice surface cause the volume to grow more slowly 
with energy; the volume scales like E^ up to 200 TeV, E^ near 1 PeV, and 
by 10 PeV the volume is increasing linearly with energy. 
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Sensitive Volume: d = 900 m: Incident Neutrino Angle = 90 deg. 
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Figure 13: Visualization of the effective volume of an antenna located 900 
meters below the ice surface. The two sets of contours outline the regions 
of SNR > 1 for shower energies of 1 and 100 PeV. At 1 PeV, the sensitive 
region already extends over a kilometer away from the antenna location. 
This volume covers different regions of ice for each possible incident neutrino 
direction. Here we show the situation for a horizontal flux incident from the 
right hand side of the hgure. 
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At first glance, the angular dependence seen in Fig.(14) is unexpected. 
Based on antenna gain, one might guess that the antenna is most sensitive to 
the horizontal neutrino flux (he. angles ~ 90°). Our results show a relative 
minimum for horizontal flux. The reason for this is the fact that the incident 
neutrino direction and that of the signal pulse differ by the 56° Cherenkov 
angle in ice, hence the maxima at roughly 90° + 9c and 90° — 9c- A more 
subtle contribution to this same effect is the favorable relative polarization 
of the signal’s electric held and the antenna height for angles signiflcantly 
above and below horizontal. 

The depth dependence can be seen in Fig.(15), where we plot the effective 
volume versus nadir angle for Eg = b PeV and four depths from 100 and 
1300 meters. Recalling that the ice is becoming increasingly radio-opaque 
with depth, and the fact that shallower deployments will be more severely 
affected by the surface cutoff, one can readily understand the qualitative 
behavior seen in this flgure. Because of Barth shadowing of the incoming 
flux, we want to optimize the sensitivity over the upper hemishpere. It is 
not clear from this flgure which depth will maximize the number of observed 
events. In the next section we will look at event rates versus antenna depth. 

The important result seen in Figs.(13) through (15) is that sensitive vol¬ 
umes of order Ikm^ are possible with as few as 1 or 2 receivers at 1 to 2 
PeV, and roughly 10 receivers at 500 TeV. These numbers indicate that un¬ 
precedented target masses for UHB detection are feasible using radio as the 
detection mechanism. This is precisely what is required if the currently pre¬ 
dicted UHB fluxes from AGN are to be observed. In the next section, we will 
fold in the expected cascade rates based on model AGN fluxes and arrive at 
deflnite rate predictions. 


10 Event Rates 

We now have in place all the elements to arrive at event rates for one of our 
proposed antennas. The effective volumes given in Section 9 can be combined 
with the BM shower rates presented in Section 5 to give a prediction for 
the number of events induced by electron and muon neutrinos. These rates 
are functions of the event energy and the incident flux angle and can be 
integrated to give a total rate. This integrated rate for the SP and ST 
AGN flux predictions is shown in Fig.(16) as a function of antenna depth. 
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Figure 15: Effective volume for observing events having an energy of 1 PeV 
as a function of incoming neutrino direction and antenna depth. Each curve 
is labeled with the receiver depth in meters. 



Observed Events Per Year 



Figure 16: Total number 
versus antenna depth for t_ 
depth is about 600 meters 




We emphasize that the rates shown are for a single antenna. Any realistic 
deployment, even in the case of a small pilot experiment, would likely employ 
an array of such sensors. The overall rate is maximized for an antenna depth 
near 600 meters. The rate per antenna is 437 events per year for the SP flux 
and 201 events per year for the ST prediction. The result using the AT flux 
as input is a completely negligible 0.06 events per year. Clearly atmosperic 
neutrinos are not a background for radio detection of AGN neutrinos. 

Figure (17) shows the logarithmic energy dependence of the event rate 
(see caption) for a hxed incident neutrino direction of 140 degrees. It is the 
product of a rising effective volume and the falling EM shower rate, also 
shown. Taking the curves corresponding to the SP prediction as an example, 
one sees that below 1 PeV, the detector volume is rising fast enough to 
‘beat’ the falling shower rate resulting in an increasing rate for observed 
events. Beyond the 1 PeV mark, the volume rises more slowly, the shower 
rate becomes more steep, and as a result, the observed rate peaks near 1 PeV 
and then begins to decrease. In the ST case, the peak occurs near 10 PeV 
because of the flatter shower spectrum. 

The energy-angle dependence of the AGN rates are shown in Figs.(18) 
and (19) for an antenna depth of 700 meters. The contours indicate the 
concentration of events coming from about 50 degrees above horizontal and 
centered near 1 PeV in the SP case and 10 PeV in the ST case. In the next 
section we discuss how the primary neutrino spectrum can be estimated from 
the location of this peak in the distribution. 


11 Spectrum Estimation 

A full scale Antarctic experiment would likely employ on the order of 100 
radio receivers arrayed along several vertical strings. Precisely how the anten¬ 
nas are arranged spatially is a detailed experimental question that will likely 
depend strongly on which measurements one wants to optimize. Assume for 
the moment that data is collected from an experiment in which the angle 
and energy of each event can be determined with enough precision to begin 
forming an image of the distribution shown in Figs.(18) or (19). According 
to the current flux estimates, many thousands of events might be collected 
yielding a reliable picture of the distribution peak as well as its rate of falloff 
in the {Eg, 9) plane. In this case, it would be a simple matter to unfold the 
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Shower Energy (PeV) 


Figure 17: The event rate (solid line) per year is the product of the effective 
volume (dash-dotted line) in km^ and the EM shower rate (dashed lines) per 
year, per km^, per steradian, per Log(energy). We show results corresponding 
to an incident neutrino direction of 140 degrees with respect to the nadir for 
the SP and ST flux predictions. The difference in spectral indices of the SP 
and ST cases results in events peaking at 1 and 10 PeV repectively. 
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Figure 18: The distribution of AGN neutrino events as detected by a 
antenna located 700 meters below the ice surface using the SP flux predi 
We plot 10 equally spaced contours of the quantity dNjdcosOjdlogE. 
the ratio between the innermost and outermost contours equal to 10. I 
are concentrated about an energy of 1 PeV and 140 degree nadir angle 
integrated rate is 437 events per year. 
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Figure 19: The distribution of AGN neutrino events as detected by a single 
antenna located 700 meters below the ice surface using the ST flux prediction. 
We plot 10 equally spaced contours of the quantity dNjdcosOjdlogE^ with 
the ratio between the innermost and outermost contours equal to 10. Events 
are concentrated about an energy of 10 PeV and 140 degree nadir angle. The 
integrated rate is 201 events per year. 
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effective detection volume, V{Es,9), from the data directly giving estimates 
of the AGN neutrino spectrum (horizontal flux, no Earth shadowing) and 
its attenuation by interactions in the Earth (angles below horizontal). The 
neutrino spectrum is critical to understanding AGN as particle accelerators. 
Important information on neutrino cross section and nucleon structure could 
come from the observed attenuation rate in the upper layers of the Earths 
crust (where target densities are known with some precision). 

One would always like to have the largest possible number of receivers. 
However valuable astrophysical information is possible even with a more 
modest array of sensors. In the case where only hundreds of events are 
measured, with energy resolution of ~ 50%, say, the distribution peak in 
energy, Emax, would still be estimated with some precision. It can then be 
used to estimate the neutrino spectrum as follows. The overall rate goes as 

R{E) oc V{E)E-^ , (14) 

with V{E) the known detection efficiency (effective volume) and 7 the un¬ 
known integral spectral index. The condition ^ = 0 applied at the distri¬ 
bution maximum yields 


dlogV 
dlogE ’ 

where we have assumed constant 7 in the region of the peak. 

In Eig.(20) we plot 7 (£'inax) using the effective volume shown in Eig.(14) 
for a nadir angle of 140 degrees, the angle where the distributions of Eigs.(18) 
and (19) peak. Eig.(19) indicates Emax ~ 8 PeV which corresponds to a 
spectrum going like 7 = 1; this is indeed the case for the ST flux at this 
energy. The SP case shown in Pig.(18) has Emax ~ 1 PeV. This corresponds 
to the spectral index 7 = 1.8 according to Pig. (20). A look at Pig. (2) reveals 
that 1 PeV is the location of a transition in the SP spectrum from 7 = 1 
below 1 PeV to 7 = 2.35 above. The 7 = 1 region ‘wants’ to have its peak 
near 10 PeV as in the ST case, and the 7 = 2.35 region ‘wants’ to peak 
at about 300 TeV. The net result is an abrupt peak at the transition point 
with our estimate 7 = 1.8 indicating an intermediate value. The information 
that the peak in the SP spectrum represents a ‘knee’ in the spectrum rather 
than a constant is contained in the shape of the distribution and cannot be 
determined from the peak location alone. 
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Figure 20: Relation between the peak event distribution energy, -Emax, and 
the integral spectral index for EM showers, 7(-E'max)- This spectral index is 
in turn closely related to the primary neutrino spectrum. 



















12 Summary of Results and Outlook 

We have presented a detailed analysis of the problem of radio detection of 
UHE cosmic ray neutrinos in Antarctic ice. This study, motivated by the ex¬ 
citing prospect of design and construction of a km^ scale detector, is the most 
thorough to date. We have included every known effect: particle physics in 
the fundamental interaction cross sections, the angular dependence of atten¬ 
uation in the Earth, electromagnetic shower evolution, coherent Cherenkov 
radio emission, polarization of the emitted held, pulse dispersion in the ice 
and antenna characteristics. Assuming primary UHE huxes from AGN as 
calculated by Stecker et al. and by Szabo and Protheroe, we calculate very 
conservatively that there will be at least 200 events recorded per detector 
element per year. Our results show the tremendous potential for the radio 
detection technique. 

Compared to earlier studies by our group and others, this work has made 
a signihcant improvement by considering buried detectors. The main advan¬ 
tage of radio detection is the long range of the radio signal in cold ice, which 
allows one to detect events efficiently more than a kilometer from the source. 
In earlier work, it was found that the Earth’s presence severely attenuates 
the interesting PeV neutrinos due to QCD effects that enhance the UHE 
neutrino-nucleon cross section. Our recent re-analysis of this cross section 
using new HERA data shows that previous calculations are more than a fac¬ 
tor two too low, so the consideration of the Earth’s screening is even more 
important than previously thought. The screening makes receivers which fo¬ 
cus strictly on upward-moving neutrinos less effective just in the region where 
radio has special advantages. By studying the response of buried detectors as 
a function of deployment depth and incident neutrino angle, we have found 
an optimum depth at roughly 600 meters where one achieves a substantial 
increase in rate. The increase in rate comes partly from including detec¬ 
tion of sidewise and downward moving primaries and partly from our better 
treatment of the geometrical factors. The background from atmospheric neu¬ 
trinos is found to be negligible, as in previous studies of downward looking, 
surface detectors. A biconical antenna design was chosen for illustration: it 
has dimensions compatible with phototube deployment and good frequency 
and angular response for our purpose, as discussed in Sections 6 and 7. A 
single antenna could detect AGN in either of the two neutrino flux models 
we use. We have also shown (Section 11) that even a modest data sample 
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with energy resolution on a logarithmic scale allows one to resolve the char¬ 
acteristic spectral index of the source as well. This technique can distinguish 
clearly between the two flux models considered here, as we have illustrated. 

We now turn to the role of muons, discussed in Section 5, and the ad¬ 
ditional effects of hadronic showers. Recall that muons contributed a small 
effect from bremsstrahlung, and that their hadronic showers were ignored in 
the calculations we presented. Our calculations, which are based on the reli¬ 
ably modeled electron-induced showers, establish a hrm lower bound on the 
capabilities of radio detection. It is not difficult to estimate the additional 
enhancement from including hadronic showers if the assumptions of Provorov 
and Zheleznykh are used. In the best case that the hadronic shower adds its 
detectable electromagnetic energy to the electron shower energy linearly, the 
electron events then approximately double their shower energy. Moreover, 
adding the muon’s hadronic shower (and still treating the muon itself with 
continuous energy loss) allows one, in effect, to add an electron shower to 
every muon neutrino induced event. Finally, including neutral currents and 
the scaling of the effective detection volume with energy, we estimate that 
each antenna would detect from 4 (Stecker et ah model) to 5 ( Protheroe 
and Szabo model) times more detected events than those we present in the 
body of the paper. To use these numbers, the reader can multiply any of 
the relevant curves by the corresponding factor: for example, the integrated 
rate for the SP flux from Fig.(18) becomes roughly 2200 events per year, per 
detector element. 

As a last illustrative exercise, we would like to present an estimate for an 
antenna array - a model for a PeV-scale neutrino telescope. There are many 
detailed considerations in array design that we do not attempt to address. 
Consulting Fig.(13), we note that approximately ten detectors arranged a 
hundred meters or so apart on a vertical string roughly a kilometer long 
would have non-overlapping volumes and linearly additive rates. This guides 
our estimate. Ten such strings separated by 1/2 km would seem to be a 
reasonable array. In a real telescope, of course, there are a number of factors 
that determine the optimum overlap between receivers. For purposes of dis¬ 
cussion, we have simply included a 2:1 redundancy on each detector volume 
- that is, a detector array with an effective volume 1/2 the value of the sum 
of the volumes of the individual, widely spaced detectors. We believe this is 
conservative, because correlations and signal processing can improve the de¬ 
tectable signal-to-noise ratio of a group of receivers to a value well below that 
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Figure 21: Summary of detected events per year using the SP and ST AGN 
flux predictions for a 100 receiver array. The volume redundancy we assume 
is 50%, that is, each observed event appears in two receivers on average. 
The dark shaded bars represent a lower bound on the observed rate based 
on pure EM showers induced by charged-current neutrino interactions. The 
lightly shaded bars show an upper bound on the rate which would occur if 
the energy in hadronic showers contributes to the signal as efficiently as the 
pure electromagnetic component. A ^itailed treatment of hadronic showers 
should yield a rate between the upper and lower bounds. 









of a single receiver, thus increasing the effective detection volume per antenna 
[35]. We do not include this effect in our estimate; in any case the reader 
can scale our results by any suitable redundancy factor in a straightforward 
way. In Fig.(21) we present the event rates calculated with such an array 
for the SP and ST models, showing both our conservative lower limits (no 
hadronic showers, few muon neutrino -induced showers) and the optimistic 
upper limit ( including hadronic showers and neutral current effects). The 
upper limits are 22,000 and 10,000 for SP and ST, respectively, representing 
good statistics that might be sufficient for creating true “images” of the PeV 
neutrino sky! 

The prospects of scientihc pay-off by exploiting radio detection of UHE 
neutrinos are excellent, and this method’s complementarity with standard 
methods is very attractive. We hope that our results will contribute to the 
push for a pilot experiment to begin the study of UHE neutrinos as soon as 
possible. 
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